Phonon-assisted excitation transfer in quantum dot molecules 
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We derive a quantum-kinetic description of phonon-assisted Forster transfer between two coupled quan- 
tum dots (a quantum dot molecule). We show that the exciton state decays to the ground state of the 
QDM via a combination of the Rabi rotation and exponential decay. For moderately spaced dots this 
process takes place on a picosecond time scale. 

Coupled quantum dots (QDs), referred to as quantum dot molecules (QDMs), have attracted 
much attention in recent years Besides the electronic coupling or superradiance effects 0, 
the properties of QDMs are affected by phonon-related phenomena, fn particular, many ex- 
periments have shown phonon-assisted excitation transfer [SJ between the QDs. For moderate 
separations between the dots 6 nm) tunneling is exponentially suppressed and the energeti- 
cally lowest states correspond to spatially direct excitons localized in individual QDs [1]. Such 
states are bound by the Coulomb interaction via interband dipole moments i-e., by the 
Forster interaction Q . Signatures of such coupling were indeed found in a photon-correlation 
experiment 

In this paper we derive a quantum-kinetic description of the evolution of an exciton in a QDM, 
including the effect of the carrier-phonon coupling. We show that, depending on the parameters, 
the system can show different dynamical scenarios, ranging from partial pure dephasing to an 
almost exponential transfer. As we shall see, the transfer may be very fast (on the time scales 
of several picoseconds), that is, 2-3 orders of magnitude faster than suggested by the existing 
perturbative estimates [S]. 

We consider two flat, cylindrically symmetric, coaxial QDs, separated by the distance D along 
the z axis and interacting with phonons. The formalism will be restricted to the subspace, 
spanned by the states |0), |1), corresponding to a single exciton in the ground state of the lower 
and upper dot, respectively (with a fixed polarization). We assume that the wavefunctions of 
excitons confined in different dots do not overlap, so that no phonon-assisted transitions are 
possible. The Hamiltonian of the system is then 

H^^Aa, + Va, + hJ2^kbibk+ J2 10(^1 E (si"^-^ + ^i'^*^) , W 

k 1=0,1 k 

where CTi are Pauli matrices in te basis (|0), |1)), A is the energy mismatch between the dots, 

V is the amplitude of the Forster coupling, 6^,6fc are creation and annihilation operators for 

the phonon mode with a wave vector k, ujk is the corresponding frequency and g^''' are exciton- 
phonon coupling constants. For heavy-hole excitons confined in QDs stacked along z one has 

V — (P/{4:TTeoei-D^) — {3/4:)[hc/{DE)]^hT, where d is the interband dipole moment, eq is the 
vacuum dielectric constant, er is the relative dielectric constant of the crystal, c is the speed of 
light and F is the inverse exciton lifetime. 

The most effective interaction between neutral excitons and phonons is the deformation 
potential coupling to longitudinal acoustic phonons. Approximating the exciton wave function 
by a product of identical Gaussians one gets Q 



where (TejCTh are the deformation potential constants for electrons and holes, v is the normal- 
ization volume for phonons, k±, are the components of the wave vector in the QD plane and 
along z, l±, Iz are the confinement sizes in these two directions, and u is the speed of sound. 
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The evolution of the interacting carrier-phonon system is described using the correlation ex- 
pansion technique (lol. [lllll^. One starts from the three dynamical variables x, y, z describing 
the carrier state, x = {ax{t)), . . ., where ai{t) = e^Ht/h^^^^iHt/h ^^^^ ^-^e operators in the Heisen- 
berg picture. From the Heisenberg equations of motion one finds the dynamical equations for 
these three variables, 

i = i([i/,a,]> = -Ay-42/^ReBfc-4y^Reyfc, (2) 

k k 

and analogous for y and z (from now on, the time dependence will not be written explicitly). 
Obviously, this set of equations is not closed, but involves the new phonon variables Bf, — gk{bk), 
as well as phonon-assisted variables of the form y^ — gk{{<7ybk)) — {f^ybk) ~ {o'y){bk)- Here 
gk = {gl^^ — gl^^)/2 and the double angular brackets, ((. . .)), denote the correlated part of a 
product of operators, obtained by substracting all possible factorizations of the product. 

Next, one writes down the equations of motion for the new variables that appeared in the 
previous step, for instance, 

ijk = i{[H,yk]) ^ Axk - 2Vzk - iujkyk + \gkfiiyz + x) (3) 

+2 ^{Xqk +Xqk)+ 4:Xk ^ReBg + 2x^{Bqk + Bgk), 
Q q q 

where the new two-phonon and two-phonon-assisted variables are defined as Bqk = 

gqgk{{hqbk)) , Bqk = g*qgk{{b\bk)) , Xqk = gqgk{{<Jxbqbk)) , Xqk = g*qgk{{(7xb\bk)) , CtC. In tllC 

next step, one writes the equation of motion for these new variables, introducing three-phonon 
variables. It is clear that the resulting hierarchy of equations in infinite and has to be trun- 
cated at a certain level. Here we do this by setting all the correlated parts of three-phonon 
and three-phonon assisted variables equal to zero. This amounts to neglecting three-phonon 
processes (that is, emission or absorption of three or more phonons within the memory time 
of the phonon reservoir, which is of order of 1 ps). The motivation for this procedure is that 
higher order correlations should play a decreasing role in the dynamics. From the equations of 
motion it is also clear that such higher order correlations develop in higher orders with respect 
to the coupling constants gk- 

In this way we treat the problem at the same level as in the recent work p^ . taking into 
account the coherent and non-equilibrium phonons which are essential for the correct description 
of carrier phonon- kitetics in QDs [l2| . As an improvement over the standard approach [ill , [l^ , 
we define collective, frequency-dependent variables 

Bl{uj) = '^5{UJ - UJk)Bk, B2{UJ,UJ') = J2qk^i^ ^ iOk)S{uj' - UJk)Bqk, 
k 

Xl{uj) ^^S{UJ - UJk)Xk, X2(UJ,UJ') = J2q^S{uj ~ UJk)S{uj' - UJk)Xqk, etC. 
k 

All the equations of motion can be rewritten in terms of these frequency-dependent variables. 
For instance, Eqs. ([2]) and ([3]) read 

X = -Ay-iy J dujReBi{uj) - 4: J dujReyi{uj), 

yi{uo) = Axi{uu) — ^Vzilu) ~ iuoyi{uo) + J{ijj){iyz + x) + Axi{uj) j dcj' ReBi{uj') 

+2 / duj'[x2{oj',u;) + X2{uj',uj)] + 2x / doj' \ B2{uj' , uj) + B2{uj' , oj)] , 
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FIG. 1: (a,c,e) The Bloch sphere representation of the evolution in the three cases desribed in the text, 
(b) The overlap of the system state with the initial state |+). (d,f) The occupation of the higher energy 
state as a function of time for two values of the Forster interaction. Here l± = 4.5 nm, Iz — 1 nm, 
D = 6 nm and the temperature T = 0. 



where J{lu) — X]fc l5fcP'^(^ ~ ^k)- In this way, the set of variables labeled by points in the 
3-dimensional reciprocal space is replaced by a set labeled by points on a real frequency axis. 
For the calculations, the frequency axis is discretized and the ordinary differential equation for 
the resulting variables is integrated numerically (we use 601 uniformly spaced points up to the 
frequency cutoff at 20 ps""'^, which yields 5.8 • 10^ variables). 

Let us begin the presentation of the results with the noninteracting case, V ^ 0. The occupa- 
tion of each dot is then conserved, so that no excitation transfer may take place. Nonetheless, 
this does not mean that no carrier-phonon kinetics takes place. In Fig. [T] we present the evo- 
lution after an instantaneous preparation of the state |-|-) = (|0) -I- |1))/V2 in a pair of QDs 
with A = 2 meV. The pure dephasing effect, related to the lattice response to the apperance 
of a charge distribution reduces the coherence of the superposition state. This is also 
manifested by the decreasing amplitude of oscillations of the projection on the initial state, 
(-|-|p(t)|-|-), where p(t) is the reduced density matrix of the carrier subsystem. 

In the presence of an interaction, V 0, the occupations of the two dots are not conserved 
and excitation transfer becomes possible. As can be seen in Figs.[Hc-e), the system evolution is 
a combination of a rotation around a tilted axis, defined by the eigenstates of the unperturbed 
Hamiltonian, and dephasing resulting from the interaction with phonons. As a result of the 
latter, at T = the system relaxes towards the lower eigenstate of the unperturbed Hamiltonian 
(with some correction due to phonon-induced energy shifts and polaron effects). For 1/ <C A, 
this lower eigenstate is close to the state In this limit, the transfer is nearly exponential, 
as can be seen in Fig. [Hd) . 

When the interaction gets stronger the final state gains more admixture of the higher-energy 
dot. Now, the transfer takes place via clearly marked oscillations about an exponential trend. 
The value of = 6 meV used in Figs. [Ue,f) corresponds to the interband diple moment of 
9 eA (or a radiative lifetime of 470 ps) in the dipole approximation, which is within the range 
of typical values for self-assembled structures. For this value, the decay rate is about 10 ps, 
which means that the process is very fast. 

The results presented in this paper show that the phonon-assisted excitation transfer in QDMs 
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is an efficient process that may considerably affect the observable properties of these systems. 

For moderately spaced dots, the Forstcr coupling results in a fast transfer showing a pattern of 
oscillations around an exponential curve. We believe that these characteristic signatures should 
make it possible to identify this process in experimental observations. 

P.M. is grateful to V. M. Axt for many discussions on the correlation expansion. 
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